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1 INTRODUCTION 



Magnetic fields are an integral part of nature. No physi- 
cal description can be complete without considering them, 
whether one deals with the realm of the very small 
(e.g. quantised Hall effect, von Klitzing 1986) or that of 
the very large (e.g. cosmological magnetic fields, Zel'Dovich 
1965). Hence, the study of magnetic fields is essential in 
understanding the universe. The birth of radio-astronomy 
in the 1930s allowed, amongst others, the measurement of 
galactic magnetic fields. The closest and most relevant to us 
galactic magnetic field is, of course, that of the Milky Way. 

The large-scale component of the Galactic magnetic 
field (hereafter GMF) is intertwined with a number of inter- 
esting astrophysical processes, such as the deflection of ultra- 
high-energy cosmic rays (UHECRs) by the GMF (e.g. Med- 
ina Tanco et al. 1998), or the existence of star-forming re- 
gions (e.g. Novak et al. 1998). A number of observables have 
been employed towards mapping the GMF: e.g. Zeeman 
splitting of spectral lines, optical starlight polarisation and 
Faraday rotation of polarised radio sources. However, the 
main challenge that most of these methods face arises from 
the observer's location inside the field's physical boundaries. 
Faraday rotation of polarised sources exploits this fact by 
providing a measure of the strength and direction of the line- 
of-sig ht (los) component of the GMF. As was first pointed 
out in iLvne fc Smith! (|l968l ) , Faraday rotation of radio pul- 
sars, in particular, is a powerful method for mapping the 
GMF. One can summarise the value of pulsars towards this 
goal as follows: (a) pulsars are distributed throughout the 
entire volume of the Galactic disc; (b) they are amongst 
the most polarised radio sources known; and, finally, (c) the 
dispersion of their pulses, given by their dispersion measure 
(DM), not only allows one to estimate their distance to a 
high accuracy, but also, when combined with their rotation 



ABSTRACT 

We measured a sample of 150 pulsar Rotation Measures (RMs) using the 20-cm re- 
ceiver of the Parkes 64-m radio telescope. 46 of the pulsars in our sample have not had 
their RM values previously published, whereas 104 pulsar RMs have been revised. We 
used a novel quadratic fitting algorithm to obtain an accurate RM from the calibrated 
polarisation profiles recorded across 256 MHz of receiver bandwidth. The new data 
are used in conjunction with previously known Dispersion Measures (DMs) and the 
NE2001 electron-density model to study models of the direction and magnitude of the 
Galactic magnetic field. 

measure (RM), provides a direct way of estimating the av- 
erage value of the interstellar magnetic field. 

The potential of Faraday rotation as a probe to the 
GMF has motivated campaigns to tr y and measu r e pul - 
sar RMs. The fi r st su ch efforts, by iManchestei] (|l972T l 
and iManchesterl (|l974T ). brought the number of deter- 
mined RMs to 38. An im portant contribution was made by 
[Hamilton fc Lvnei (|l987l ). who measured 163 values, increas- 
ing the total nu mber of avai l able RMs to 185. Sub s equen t 
observations by Costa et all dl99ll). iRand fc Lvnel 119941) 
lOiao et all (|l995r T Ivan Ommen et all (|l997l ) and lHan et al.l 
( 1999) brought the total to 320. To date, the pulsar cata- 
logue contains 554 RMs: i.e. nearly a third of the pulsars 
known have an RM value. 

The first use of available pulsar polarisation data to 
map the GMF led to unequivocal evidence for a clock- 
wise (CW) directed local field, as viewed from the Galac- 
tic north, of uniform strength ~ 2 fiG (Manchester 1972; 
Manchester 1974). Re - analy sis of 48 pulsar RMs from the 
[Manchester fc Tavlorl (|l977l ) catalogue provided the first 
indication for a magnetic-field reversal in the Carina- 
Sagittarius arm (jThomson fc N elson 198Q ) . Thes e finds were 
confirmed by the work of ILvne fc Smith! ([1989 ), who also 
suggested an additional reversal outside the Perseus arm. 
As more data were being collected, t he field's structure to 
larg er distances could be explored: e.g. IRand fc Lvnd (|l994l ) 
and lHan et ail l|2002l ) reported magnetic-field reversals in 
the Crux-Scutum and Norma regions. 

Despite the wealth of information that observational 
data can provide, the sole use of pulsar RMs and DMs for 
the estimation of the interstelar field cannot provide suffi- 
cient information about the GMF structure: this is because 
the quantity RM/DM can only provide a weighted average of 
the parallel component of the field, with the averaging per- 
formed across the distance between pulsars. But since both 
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the electron density and the magnetic field are far from sim- 
ple functions of azimuth and distance, inverting RM/DM to 
obtain B is not feasible. 

A possible solution is to try and fit magnetic-field mod- 
els to the data and find which model is more representative 
of the observed RMs (see section 6). A number of authors 
have tried vari ous field configura t ions to match the avail- 
able data: e.g. ISofue fc Fuiimotol (|l983h were amongst the 
first to favour a 2-arm bi symmetric-spiral mo del, whereas 
IRand fc Kulkarnil |l989l) . iRand fc Lvnel (|l994l ) and IValleei 
I 2005h iitted a concentric ring model: the latter type of 
model showed consi stency with earli er data, but more re- 
cent observations bv lHan et al.l (|2006h contradict the model 
predictions, especially towards the inner Galaxy. 

Despite those efforts, however, the shape of the GMF re- 
mains unclear. The main difficulties in modelling arise from 
(a) a short-scale, seemingly random magnetic-field compo- 
nent — of the same order of magnitude as the large-scale one 
— that causes large RM fluctuations across the entire Galac- 
tic plane (hereafter GP); (b) regional anomalies of the ISM 
(e.g. the Gum Nebula, the North Polar Spur) , which corrupt 
a potentially smooth RM variation; and, in addition, (c) the 
pulsar distances — upon which a large portion of modelling 
is relied — are not certain enough, which causes discontinu- 
ities in the field direction over large longitude ranges. 

In this and forthcoming publications we will study the 
difficulties of modelling the GMF in some detail. In this 
paper we will concentrate first on the measurement of new, 
accurate RMs in order to supplement previously sparsely 
sampled regions of the Galaxy. We will then examine as to 
whether commonly studied models of GMF are consistent 
with the new constraints provided by our often more precise 
measurements, and how they can be improved to provide a 
better fit to the data. A full modelling of the GMF with the 
full sample of available data will be undertaken in a later 
publication. 



3 DATA ANALYSIS 

We proc essed all the data using the PSRCHIVE software 
package (|Hotan et al.ll2004h . within which we developed new 
software to perform RM fitting. The raw data from the wide- 
band correlator first underwent a process of RFI excision and 
polarisation calibration. Following this, each data set was 
summed in time to obtain four integrated Stokes profiles for 
each frequency channel. One average polarisation position 
angle (PA) was computed per frequency channel from the 
Stokes parameters Q and U, summed across n phase bins i 
corresponding to the pulse, as: 



PA = o arctan I v-n n 



(1) 



We computed the error of this average PA using 
one of the two following methods, depending on the 
signal-to-noise ratio of the total linear polarization L = 
(IXo^ + E^ 2 ) 172 , given by s/n = (l/^)L/ar, 
where 07 is the RMS of the total intensity. For low 
s/n ratios (< 10), the la error on the PAs, <jpa, was 
computed numerically using the method described in 
iNaghizadeh-Khouei fc Clarke! (|l993l ) . For high s/n PAs, i.e. 
^ 10, the PA error was c omputed using the formula from 
lEverett fc Weisberd (|200lD : 



apA = ~T 



(2) 



In obtaining the average PA across the pulse per fre- 
quency channel as described above, there is a possibility that 
the integration of Q and U may in fact lead to a low average 
linear polarization and, therefore, a large error on the PA. 
This depends on the shape of the PA profile; integration of 
consecutive bins with orthogonal PAs significantly reduces 
the total linear polarization. While this possibility exists, we 
found that in our data the error on the average PA was sig- 
nificantly improved with respect to the errors on the PA in 
individual bins. These reduced errors were instrumental to 
good RM fits, therefore we deemed the process of averaging 
the PA appropriate. 



2 OBSERVATIONS 

We observed a sample of 239 pulsars with the 20-cm H-OH 
receiver of the Parkes 64-m radio telescope in order to obtain 
a polarization database at high time resolution for all these 
objects. The observing session was carried out from 2006 
August 24 to 27. The H-OH receiver at Parkes is equipped 
with two, orthogonal linear feeds that are sensitive to po- 
larised emission with an equivalent system flux density of 
43 Jy (Johnston 2002). The full Stokes profiles of each pul- 
sar were recorded using digital correlator back end, across 
256 MHz of bandwidth split in 1024 frequency channels. The 
wide-band correlator is capable of producing high-resolution 
profiles with 1024 phase bins, and more than 60% of the pro- 
files were recorded in this mode, while the rest had 512-bin 
resolution. We observed each of the 239 pulsars for ~15 min, 
producing ~90 10-s subintegrations. In order to save disk 
space, we further time-averaged the data to reduce them 
to 15 1-min subintegrations, which we later calibrated and 
analysed. 



3.1 The Rotation Measure from a Quadratic Fit 
with Wraps 

The degree of Faraday rotation, APA, that electromagnetic 
waves of wavelength A undergo as they propagate through 
the interstellar medium from the pulsar to the Earth is ex- 
pressed by the Rotation Measure (RM). This rotation is 
caused by the interaction of the radio waves with the mag- 
netised plasma of the interstellar medium and, more specifi- 
cally, the plasma electrons along the line-of-sight (los) to the 
pulsar. If n e is the number density of the plasma electrons 
along a unit column of length equal to the distance to the 
pulsar, d, and B, the magnetic field of that plasma, then 
RM is given as the path integral of n e B along the line of 
sight: 

RM =^f ^ S)B[S) - dS (3) 
where ds is the path vector element in the direction of wave 
propagation, e is the electron charge, m e , the electron mass 
and c, the speed of light. Using cgs units for the constants 
and expressing the distance in pc and B in fiG, we get 
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Figure 1. Examples of RM results obtained using a quadratic fit 
with wraps. The position angle, PA, is plotted against frequency. 
The solid lines represent the best fit from the method described 
here, whereas the dashed-dotted lines corres pond to the fits based 
on previously published RMs bv lHan et al. (2006). 



RM = 0.812 



n e (s) 



B[s) 



fiG 



ds\ 
pc) 



(4) 



so that RM is expressed in the usual units of rad m -2 . 

An observable property of the above magneto-optical 
interaction is that, for a particular pulsar, the longer the 
wave, the greater the effect of Faraday rotation, with in- 
finitely short waves being unaffected. Mathematically, this 
realtionship is expressed as 



3.2 Statistical Error on the RM from Monte 
Carlo Simulations 

In order to estimate the statistical error on the best-fit RMs, 
ffstat , we generated a look-up table of RM errors correspond- 
ing to the typical (average) PA error in our data, OpA ta , and 
the actual RM value that came out of the fit to the data, 

RM data 

In order to create the look-up table, we generated 
Monte Carlo (MC) data sets that contained fake PA val- 
ues which were calculated from a grid of chosen RM val- 
ues, (RM mc ), using Eq. [S] The observation bandwidth 
and central frequency in our calculations were matched to 
those of the real data. The chosen grid of RM values was 
(RM mc ) = {0, 5, 10, 20, 50, 100, 200, 300, 500, 800, 1000} rad 
m -2 . In each Monte Carlo data set, the fake PAs were the 
products of randomisation of the calculated PA (from Eq.[S} 
within the window [— OpA, "pa]- We ranged the value of a'p^ 
between 1 and 63 degrees, in 1-deg steps. For each combina- 
tion of o-pA and (RM mc ) values, we generated 2,000 Monte 
Carlo data sets. Each of the generated MC data sets were 
then fitted using the fitting routine that was used with the 
real data. The resulting RM values from the fits, RM mc , were 
Gaussian distibuted around the value (RM mc ) that was used 
to generate that data set. For each distribution, we kept the 
RMS spread as the (J sta t corresponding to the particular pair 
of {(RM mc ), o"pa} values that was used to generate the dis- 
tribution. 

The look-up table that was generated in this way al- 
lowed us to assign a statistical error to the RM data values of 
the real data, given the typical er PA ta of the PAs, by matching 
- as closely as possible — the values of cr PA ta and RM data 
to those of (TpA and (RM mc ). For values lying between the 
tabulated ones, we used a two-dimensional cubic-spline in- 
terpolation of the table values. 



APA = RMA 



P 



(5) 



where / is the frequency of the waves. 

Eq. [5] implies that the RM can be measured across the 
observation band by fitting a quadratic function to the PAs 
of every frequency channel. Using a Bayesian approach, we 
developed a fitting algorithm that was applied to the avail- 
able data sets. One of the main strengths of the quadratic 
fitting algorithm is that it accounts for the 180° ambigu- 
ity that the PAs are subject to according to Eq. [1] The 
algorithm finds the best fit to the data by means of a 
Bayesian likelihood test. The fitting function was of the 
form, PA = PA + c 2 RM(l// 2 - l// 2 ), where PA is the 
PA at frequency /o, and fj is the frequency of channel j. 
The free parameters PAo and RM were stepped through the 
ranges [0, ir] rad and [—1000, 1000] rad m -2 , in 1° and 1-rad 
m -2 steps, respectively. 

The accuracy of the method is clearly displayed in the 
examples of Fig. [1] Our confidence in the calculated RMs 
is drawn both from the well calibrated data and the robust 
fitting approach, but also from the meticulous determination 
of the RM errors. The latter, we believe, has produced more 
firm results compared to previous efforts. The following two 
sections describe the procedure for the calculation of these 
errors. 



3.3 Systematic Error Determination 

We made an effort to try and determine the systematic ef- 
fects on the data, which can be caused e.g. by imperfect po- 
larisation calibration or RFI mitigation. The procedure fol- 
lowed was that above, but at the start of each fitting process 
we averaged the number of frequency channels by a factor 
Kf, whose value we varied before each fit. More specifically, 
three different factors were used, 32, 64 and 128, which guar- 
anteed a minimum of 8 frequency channels in the fit — thus 
avoiding low statistics. It is expected that increasing the 
value of Kf will lead to increased s/n per frequency chan- 
nel; but it also means that there are fewer degrees of free- 
dom in the fit. Thus, the resulting value of the RM should 
stay the same, unless there are systematic effects that al- 
ter each result. The magnitude of those systematic effects 
can be estimated as the bias-corrected standard deviation 
of those three values. Prior to this estimate, each of the val- 
ues was weighted by the inverse square of the goodness-of-fit 
probability, pa t , arising from the Bayesian fitting. Hence, we 
calculated the systematic error as 



- (RMk / -(RMk / ))' 



(6) 
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Figure 2. Residuals' plot between previously published RMs 
(RM pub ) and those calculated from our analysis (RM ncw ), for 
108 pulsars, sorted by increasing Galactic longitude (£). The or- 
der number for each pulsar is shown in parentheses, in column 1 
of Table [T] The RM residuals were calculated by subtracting the 
published from the new RM value and dividing the result by the 
standard deviation of the new RM. The two arrows at the top of 
the plot correspond to RM measurements for PSRs J1633— 5015 
(leftmost) and J1707— 4053 (rightmost), for which the residuals 
were outside the plot's range. 



where ^RMk ; ) is the weighted mean with weights equal to 
WK f = 1/pif 



4 RESULTS 

Of the 239 pulsars in our sample, our analysis produced RM 
values for 150. The remaining 89 pulsars had either low s/n 
generally or very little linear polarization. Of the 150 RMs 
obtained, 46 are values measured for the first time and 104 
revised from their previous values. The sample of previously 
unpublished RMs covers over 200° Galactic longitude and is 
mostly confined within a 4 kpc radius from the Sun. Apart 
from being a valuable addition to previous RM databases, 
the new sample enriches our knowledge of RMs in the first 
quadrant in particular (0° ^ £ < 70°), which was previously 
poorly sampled. All resulting RM values and the respective 
combined errors (ctrm = (cr| ta t + ^sya) 1 ^ 2 ) are listed in Table 

m 

Fig. [2] shows a comparison between our sample of RMs 
and previous results from the literature. The results agree 
within the error bars for the majority of the sample. Some 
objects show disagreement at the level of up to 10 rad m~ 2 . 
This could simply be due to the measurement techniques 
used and/or an incomplete estimate of the error bars. How- 
ever, it is known that systematic changes in RM at this sort 
of level can occur over periods of years as the pulsar tra- 
verses through the interstellar medium (e.g. in the Vela pul- 
sar, Johnston et al. 2005). These differences would be worth 
following up in more detail. A few pulsars have extremely 
discrepant RMs, the nature of which is not entirely clear 
but may be as simple as typographical errors in the orig- 
inal reports. In particular, the RMs of PSRs J0255— 5304, 
J1633-5015 and J1707-4053 are wildly at odds with liter- 
ature values (see Fig. [2}. 



5 MAGNETIC FIELD MODELLING 

The value of the rotation measure for each pulsar from ob- 
servations can be used to provide an estimate of the average 
magnitude of the component of the regular Galactic mag- 
netic field that is directed along the los to the pulsar: 



Jg n e (s)B^ds 
Jo Ms) 



(7) 



where the denominator is the column-density integral to the 
pulsar, called the dispersion measure, DM. Hence, 



= 1.232 



RM 



rad rrr 



DM 



pc cm~ 



(8) 



resulting in a value for Bu in /xG. 

The form of the Galactic free-electron density, n e (s), 
is derived from modelling. In this work, we assumed the 
density distribu t ion d escribed in the NE2001 model by 
ICordes fc Laziol jl200ll) - other model s include e.g. the 
TC93 model by iTavlor fc Corded (|l993h . Furthermore, in 
what follows, we will use the Galactocentric Cartesian co- 
ordinate system of the NE2001 model, shown in Fig. 1101 
All radii in our modelling correspond to Galactocentric dis- 
tances, with the azimuth angle, 9, measured clockwise from 
the Galactic centre-Sun direction. 

For the observed sample of pulsars, the value of \Bu) is 
shown in column 6 of Table [T] Despite the fact that pulsar 
RMs allow us to probe only the average magnitude of the 
GMF along the los, the numerous, scattered pulsar positions 
across the Galactic disc — as well as a number of them in 
the Galactic halo — makes up for a powerful method to test 
the various theoretical models that describe its shape and 
magnitude. 



5.1 Regular GMF models 

Below are four analytical models which we chose to test 
for consist ency against the data. The first three are de- 
scribed in iKachelriefi et all {2007) and have the common 
component of a logarithmic-spiral disc field (first proposed 
by Simard-Normandin & Kronberg 1980). The last model 
is the pure dipolar-toroidal field (produced by an AO dy- 
namo), which is considered the most likely configuration for 
the magnetic field of the halo (e.g. Han et al. 1997, Han 
2002). In the description of the spiral models, we have indi- 
cated in the subsection title whether the model has a bisym- 
metric spiral (BSS) or antisymmetric spiral (ASS) configu- 
ration and whether the field is antisymmetric (indicated by 
"-A") or symmetric (indicated by "-S") with respect to the 
GP (z = 0). 

Previous efforts to check the c onsistency betwee n the 
available models and the data include lHan et al.1 (|l999f ). who 
concluded that an ASS with no reversals fits well pulsars at 
high latitudes (j b j > 8°), whereas the disc field can be better 
described using a bisy mmetric-spiral configura tion (BSS); a 
concentric-ring model dRand fc Kulkarnll 19891) was th e least 
supported by the data. Ilndrani fc Deshpandd |l999) have 
also tried to fit 131 pulsars to a concentric-ring and a BSS 
model: the latter gave the best overall fit. Nevertheless, it 
is worth mentioning that both models supported a stronger 
field in the interarm regions compared to that in the arms 
and that the general conclusion was that the field is reversed 
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Table 1. The results from the herein analysis (column 2), including the combined statistical and systematic errors (column 3). The 
table is sorted by the RA of each pulsar. For compatibility with Fig. [2] the ascending order according to I, of pulsars with previously 
published RM values, is given next to each pulsar's J2000 name (Column 1). Columns 4 and 5 show the RM values and their error that 
were previously published for these pulsars; the relevant publications are letter-coded, with the corresponding reference listed below the 
table. The last two columns, column 6 and 7, list the average value of the parallel component of the magnetic field, and its error, based 
on Eq. \E\ 
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as one crosses the boundary between the arm and interarm 
regions. 

5.1.1 TT Model (BSS-A) 

T he first model we explored was based on the model 
of iTinvakov fc Tkache^l (|2002t ) (TT). TT developed their 



model primarily to explain the trajectories of UHECRs from 
Active Galactic Nuclei (AGN) deflected by the GMF. Its 
shape resembles the spiral shape of the distribution of the 
optical matter in the Milky Way. The magnetic field in cylin- 
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drical coordinates is 



b(r) cos 



In 

tanp \ro 



/(*)* 
x (sinpe r + cos peg) 



(9) 



where p is the pitch angle, which describes how tightly 
wound the spiral arms are: the published value for this model 
is p — —8°. The radius ro represents the Galactocentric dis- 
tance along the line between the Sun and the Galactic centre 
(hereafter GC) to the field maximum, | B(r, 0) |. TT ex- 
press ro in terms of the heliocentric distance to the closest 
magnetic-field reversal, do: ro = (r© + do) exp(— tanp), 
where do = —0.5 kpc. A negative do implies that the clos- 
est reversal occurs between the Sun and the GC (t — 0°), 
whereas positive values of do place the closest reversal to- 
wards the Galactic anticentre (I — 180°). The functions 
b(r) and f(z) are characteristic of such "spiral" magnetic- 
field models and respresent the radial and vertical profile of 
the field's magnitude, respectively: for the TT model, 6(r > 
r CO vc) = b(r®) x (r /r), where &(r©) = 1.4 /iG, and b(r < 
7"core) — const, where r coro = 4 kpc. The vertical decay is 
given by the exponential law f(z) = sign(z) exp (— | z | /zo), 
where it is typically assumed that the value for the scale 
height of the Galactic halo, Zo, is 1.5 kpc. 



5.1.2 HMR Model (BSS-S ) 

Another spiral model w hich we tested again st our data was 
based on the model of lHarari et~aT1 l|l999l ) (HMR). This 
model shares the same magnetic-field shape in the Galac- 
tic disc with that of the TT model but differs in the field 
symmetry with respect to the GP, the latter being symmet- 
ric for the HMR model. 

HMR describe their spiral as somewhat less tightly 
wound than that of TT, with p = —10°, whereas the first 
magnetic-field reversal happens at do ~ —0.5023 kpc. The 
radial magnetic-field decay in this model is described by 
the hyperbolic-tangent law b(r) — (3r©/r) tanh 3 (r/r coro ), 
where r coro = 2 kpc. The vertical suppression in the HMR 
model follows the two-scale law /(z) = l/[2 cosh(z/zi)] + 
l/[2 cosh(z/z2)], where the scale heights z\ and Z2 are typi- 
cally set equal to 0.3 and 4 kpc, respectively. 



where -Bt, max (r©) = 1.5 /iG is the field maximum above the 
GP; fir = 1-5 kpc is the height of the field maximum and 
Wt = 0.3 kpc is the half-width of the Lorentzian field profile; 
finally, H(x) is the Heaviside step function. 

A final, dipolar component is added to the above two, 
with Cartesian components 

B x = — 3^- cos 4> sm 4> sm & 
R A 

B v — — 3^7 cos d> sin <b cos 6 
R 3 

B z = g(l-3cos 2 <« 

(11) 

where R = (r 2 + z 2 ) 1/2 , 9 = tan _1 (x/j/) and <j> = 
cos -1 (z/J?) are the Galactocentric spherical coordinates; 
/iG = 123 /j,G kpc 3 is the magnetic dipole moment. The 
singularity at R — is avoided by setting B z = —100 /iG 
for R < 500 pc. 



5.1.4 Dipolar-Toroidal Model 

In addition to the spiral models, we used a pure Dipolar- 
Toroidal magnetic field (see e.g. Han et al. 1997). The total 
field, which is the sum of the toroidal and the dipolar fields, 
is given in spherical coordinates as 



B = -Bdipole + Bt 



toroidal — 



— ^3 (2cos#e_R + sin 9eg) 



■ ^ n ~ (12) 
sign(z)— e e 



The parameters m and n represent the relative magnitude 
of the dipolar and toroidal field components, respectively, at 
a given point in the Galaxy, weighted by the Galactocentric 
distance of that point. Recent work has reported a value of 
m = 245^56 pG kpc 3 and n = 4.8lo;® A« G k P c via X 2 fitting 
of pul sar RMs from a sa mple of nearby pulsars with d < 0.5 
kpc (| Athanasiadisll2004h . 

Apart from the dipolar field of the PS model, none of 
the above model fields includes a vertical component. This is 
characteristic of most GMF models, although it is believed 
that a dipolar structure a t the GC may be t he reason for a 
weak B z w 0.2 ^G at r l|Han fc Qiaolll994l '). 



5.1.3 PS Model (BSS-S) 

The last spiral model that we used is a variation of the 
BSS-S model by iProuza fc Smfdal (|2003l ) (PS). The shape 
of the disc is identical to that described in the TT model, 
the only differences being the smaller halo size, with scale 
height zo = 0.2 kpc, and the higher field magnitude, which 
we normalised to 6(r©) = 2 /iG. 

In addition to the disc and halo components, PS con- 
sidered a toroidal contribution, Bt, in the form of circular- 
discs either side of the GP, with the vertical profile of the 
field following a Lorentzian distribution: 

B T = -sign(z)_BT,max(r0)x 



x g(r s -r) + g (r- r ? )e-o (cog ^ - sin ^) ( 10 ) 



-h T \ 

"T j 



6 THE LARGE-SCALE STRUCTURE OF THE 
GALACTIC MAGNETIC FIELD REVISITED 

Despite the conclusive evidence for the large-scale GMF di- 
rection within d ~ 1 kpc (e.g. the CW direction of the lo- 
cal field and the field reversal in the first quadrant (Ql) 
of the Sagittarius-Carina arm, Lyne & Smith 1989), past 
analyses of the available database of pulsar RMs have re- 
sulted in contradicting opini ons with respect t o the field's 
structure at larger distances: iMitra et al.l 1120031 ) studied 45 
pulsar RMs towards the Perseus arm and concluded that 
no magnetic field reversal is evident to a 5-6 kpc distance; 
more recently, a broader analysis of 223 southern pulsars 
by Han et al. concluded that the arm regions maintain a 
counter-clockwise (CCW) magnetic field whereas the data 
suggest that the field is clockwise (CW) in the interarm 
regions (|Han et alj |2006f). In addition to the usefulness of 
Galactic pulsars in studying the GMF, a number of au- 
thors have included extragalactic sources in their analyses. 
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Figure 3. (a) Hammer-Aitoff (edge-on) and (b) Galactic-plane (face-on) projections of the observed sample of pulsars which were 
used in our analysis. The Galactic-plane projection also includes a gray-scale map of the free-electron density distribution according 
to the NE2001 model. The associated symbols correspond to 150 good-quality RM measurements. Each symbol's linear dimension is 
proportional to (| RM I) 1 / 2 derived from each measurement: positive RM measurements are represented with squares; negative ones, 
with circles. The largest square corresponds to RM = 848.3 rad m -2 and the largest circle, to RM = —605.9 rad m — 2 . In addition, the 
plot further classifies the total sample into RMs from pulsars without a previous measurement (green), values from pulsars that had an 
RM recently published by Han et al. (blue), and values from pulsars with previously published RMs by other authors (red). 



In 2003, iBrown et aD (|2003l ) combined 262 RMs of extra- 
galactic radio sources from the Canadian Galactic Plane 
Survey (CGPS; Taylor et a l 1999 ) with 12 RMs of pul- 
sars observed bv lMitra et ail (I2003D and others, to conclude 
that the re are no reversals beyond the solar radius. Very 
recently, IBrown et~aH (|2007l ) measured the RMs of 148 ex- 
tragalactic radio sources from the Southern Galactic Plane 
Survey (SGPS; McClure-Griffiths et al. 2 005) and combine d 
them with 120 k now n pulsar RMs from lHan et al.1 {2006), 
lHan et al l (|l999l ) and lTavlor et~afl (|200Ch . within the SGPS 
region (253° ^ I ^ 358°). They concluded that no reversals 
exist beyond the radius of the Sagittarius-Carina arm, in 
the region surveyed by the SGPS, and that the field main- 
tains a CW direction in that region. However, they found 
clear evidence of a reversal in the interarm region between 
the Sagittarius-Carina and Scutum-Crux arms, where the 
field is CCW. 

In order to shed more light on the problem of field rever- 
sals beyond the solar neigbourhood, we analysed the sam- 
ple of 150 pulsar RMs shown in Table [1] Amongst them, 
there are 57 previously unpublished values that correspond 
to 52 pulsars, of which more than 90% are distributed within 
d = 4 kpc, but with a longitudinal spread covering more that 
200° across the sky. The longitudinal and latitudinal pulsar 
distribution in the sky is graphically presented in Fig. [3^,, 
which shows a Hammer- Aitoff projection of all pulsar po- 
sitions and their relative RM magnitude. Furthermore, a 
two-dimensional projection of the GP in Fig. [3p shows the 
radial distribution of the pulsars. 

A straightforward method to assess the radial profile 
of (-B||) is to plot its value against distance, alon g different 
lines-of-sight. According to lLvne fc Smith! (|l989l ), the mag- 
netic field in the region between distance di and d-z along 
the los can be assessed from the ratio of the RM and DM 
gradients over that region: i.e. 

ARM 



= 1 232- 
dX-d2 ■ ° ADM 



(13) 



distance, d. In practice, this method is limited by the number 
of pulsars lying within a solid angle of a certain longitude 
and latitude width, Al and A6, respectively. However, since 
most pulsars in our sample lie within | z \= 2 kpc of the GP 
(see Fig. [4^,), and because, according to the B-field models 
described above, the disc field does not extend beyond that 
z-height, it was more appropriate to plot the (By) function 
inside the zone defined by | z ^ 2 kpc, instead of using a 
latitude constraint. 

In principle, such plots can be drawn for a number of 
longitude sectors across the sky. One convenient way, how- 
ever, to map the locations of the field reversals in a single 
plot is an image map of (-By) which is constructed by pro- 
jecting the pulsar locations onto the Galacting plane and 
binning the average field values, along with the pulsar co- 
ordinates, in a polar grid (referenced at the Sun). A linear 
interpolation of the available field values is also introduced 
to smooth the GMF gradients between bins. This projec- 
tion introduces an uncertainty in the z direction, which also 
propagates into the calculation of (-By) through the scale 
factor 7 = 1/(1 + (z/p) 2 ) 1/2 , where p = (d 2 - z 2 ) 1/2 is the 
polar radius. So, for a planar GMF, B(x,y), the projected 
value of the average field expressed by Eq. [7] is modified as 
follows as one moves away from the GP: 



<*!> = 



JJ d n e (s)ds 



(14) 



Ideally, one would need an infinite number of pulsars 
in order to represent the actual fluctuations of (-By) d with 



where it is also assumed that the field direction does not 
change between the z-boundaries. Clearly, the projected 
field values are representative of the actual ones as long as 
z is chosen small enough for 7 w 1. 



6.1 Location of Reversals 

Fig. [5] shows image maps, as described above, of the pro- 
jected value of the parallel magnetic-field component binned 
along 72, 5° longitude sectors and 1,000, 10-pc distance bins. 
The number of pulsars per longitude sector is distributed as 
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Figure 4. The distributions of (a) the ^-height and (b) the Galac- 
tic longitude for the pulsar sample used in the herein analysis. The 
light-gray longitude bins correspond to the four 5° sectors that 
were checked for field reversals (see section \6.l\i . 



changes were probed to only as good a resolution as the data 
spacing allowed. 

We can comp are the above statements with those of 
iBrown et al. (2007), who provided convincing evidence for a 
CCW field inside the arm region, along £ w 307° . In general, 
our conclusions support the existence of a field with CCW 
direction inside the Crux arm. If we exclude the positive 
value at as 8 kpc (« 3 /iG), which produces a sudden reversal 
in the field, then the two analyses are fully compatible with 
a CCW field inside the arm. In addition, the presence of a 
reversal in the interarm region is drawn from both analyses. 

Finally, we used the analytical GMF functions of section 
!5.1l to calculate the discrete values of (B||) d at each pulsar's 
position and compare it with the corresponding data value. 
The results are plotted with different symbols for each model 
in Fig. [6] In addition, the continuous profiles of the models 
for 2 = show the fluctuation of the field along the respec- 
tive sector. In general, the model values are not representa- 
tive of the data and seem incapable of describing the rapid 
fluctuations of the field. This is perhaps an indication that 
the regular field component can only partly account for the 
observed field, whereas an additional turbulent, small-scale 
field is needed for the full description. For pulsars below the 
CP (Fig. |6]3) with d < 2 and > 5 kpc, one can note the 
relatively good description of the data by the PS model. 



shown in Fig.[4]D. Only bins with available measured or inter- 
polated data are shown: the interpolation was only possible 
if a bin lay between measured values; extrapolation was not 
performed. The different pixel intensities correspond to dif- 
ferent values of the average field: the graduated shades of 
blue correspond to the range of positive GMF, i.e. directed 
towards the Sun, and the graduated shades of red corre- 
spond to negative GMF, i.e. directed away from the Sun. 
Both ends of the colour scale represent a minimum absolute 
value of 3 [iG — meaning that stronger fields are represented 
by the same colour. In order to test the hypothesis of mod- 
els supporting a field reversal as one crosses the GP (i.e. the 
class of models denoted with the suffix '-A'), we plotted the 
image maps based on (a) all available data in each sector, 
(b) only pulsars that lay below the GP and (c) only pulsars 
that lay above the GP. 



6.1.1 Crux 

One of the most evident reversals revealed by the map of 
Fig. [5^ coincides with the los tangential to the Crux arm, 
covering 305°-310° of Galactic longitude. The scatter plot 
for that sector presents a complex picture (see Fig|6ji): the 
field is reversed from CW to CCW somewhere between the 
Carina and Crux arms (2-4 kpc); an increase of the field in 
the Crux arm region (5-10 kpc) is observed at larger dis- 
tances, which however fluctuates around zero on a distance- 
scale of ~ 2 kpc. The field behaviour becomes less complex 
if one separately examines the field profiles below and above 
the GP (see Fig. |B]d,c): in both z > and z < plots, the 
field direction between the Carina and Crux arms is CW; 
deep inside the Crux arm region, at ft; 7 kpc distance, the 
plots for z > reveal a second field reversal, that corre- 
sponds to a change in the field direction from CCW to CW. 
One has to bear in mind, of course, that the field-direction 



6.1.2 Carina 

The sector tangential to the Carina arm (I — 280°-285°) is 
plotted in Fig.|BJl,e,f. When all pulsars in that sector are con- 
sidered, the field dips near w 4 kpc to become CCW, where 
elswhere it maintains a CW direction. Based on the only 
two pulsars below the GP, the field magnitude is decreased 
between as 1 and 8 kpc, while maintaining a positive value. 
This profile is consistent with a field pattern that follows 
the Carina arm in a CW direction. Surprisingly, the picture 
is reversed for the pair of pulsars lying above the GP. The 
RM and DM gradients between d as 4 and 9 kpc imply a 
field whose magnitude increases from a very weak (~ —0.2 
/xG) negative value, near 4 kpc, to a much stronger (~ 1.7 
fiG) positive value, near 9 kpc, well beyond the Carina arm. 
This profile is indicative of a CCW field reversi ng to CW 
somew here between 4 and 9 kpc. In comparison, lHan et al] 
(2006) reported a CCW field in the Carina arm, but noted 
an evident change in the slope of ARM / ADM between 3 and 
5 kpc, which they referred to as the "Carina anomaly" . Such 
a reversal, which the authors support is possibly attributed 
to H11 regions, is consistent with the profile of Fig. |BJl. 

Apart perhaps from the most remote pulsar in this sec- 
tor, the comparison between models and data in Fig. [SJi 
seems to favour the Dipolar-Toroidal model as the best de- 
scription of the field fluctuations. Nevertheless, the model 
values are still more than 1 standard deviation away from 
the data, and furthermore the individual plots for either side 
of the GP show a clearly inconsistent picture between the 
models and the data points. A mitigating factor in this case 
could be the existence of Hn regions, which — if real — pose 
a clear restriction to large-scale modelling. 
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Figure 5. Image maps of the projected (By) values for different longitude sectors on the GP. The field values were binned in 72, 
5° longitude sectors and 1,000, 10-pc distance bins. Prior to calculating each pixel's value, based on Eq. 1131 we performed a linear 
interpolation of the available data. Hence most pixels do not correspond to any (projected) pulsar's position. Only bins containing a 
value based on the above interpolation scheme are plotted, and thus each sector terminates at the location of the last pulsar available in 
the sector. The pixel colour and intensity along the radial paths (each representing a longitude sector of SI = 5°) correspond to different 
(S||) d values and are graduated from red, for fields directed away from the observer, to blue, for fields directed towards the observer, 
according to the scale shown below the plots. The colour scale has been intentionally saturated at both ends to reveal the field gradients 
in —3 $C (^||) ^ 3 fiG; hence, the colour at both edges corresponds to an absolute field magnitude of at least 3 /xG. In order to reveal 
potential differences in the GMF either side of the GP, we produced the maps for (a) all pulsars, (b) those with z < and (c) those with 
z > 0. 
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Figure 6. Scatter plots of the projected (By) as is calculated from the gradients of RM and DM between the available pulsars lying 
along three longitude sectors: 305° ^ i ^ 310° (top row), 280° < £ < 285° (middle row) and 30° < £ ^ 35° (bottom row). For each 
longitude sector, the field profile is shown based on the entire sample of pulsars (left column), only those that lie below the GP (middle 
column), and only those that lie above the GP (right column). The projection was applied to the sample of pulsars inside — 2 sC z 2 
kpc. A magnetic-field reversal in the region with 3 SC d sC 5 kpc, coinciding with the far edge of the Carina-Crux interarm region, is the 
most likely explanation for the sudden transition from positive (field directed towards the Sun) to negative values of the projected field. 
Also plotted are the model values of (By ^ from each model at each pulsar's position: red triangles for the TT model, green circles for 
the HMR model, blue stars for the PS model and purple open circles for the Dipolar— Toroidal model. Finally, the continuous function 
of <^_B|| ) rf with distance, for each model, is plotted for 2 = 0: the solid lines correspond to the central longitude value of i = 307.5°, and 
the dashed lines, to the edges of the longitude region. 
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6.1.3 Scutum 

Another interesting longitude sector on the image map of 
Fig. [5^ is the one that is tangential to the Scutum arm, 
between 30° and 35°. It is also a direction which was pre- 
viously poorly sampled and benefited mostly from the new 
RMs. According to Fig. [6g, a reversal occurs between 1 and 
2 kpc distance, inside the Sagittarius arm. The reversal is 
located near the far edge (from the Sun) of the Sagittarius 
arm, and it changes the field direction from CW to CCW. 
Fig. |6j>',h seem to support a predominantly CW field in the 
Ql part of the Sagittarius arm, although more pulsar RMs 
are needed in that region to have a solid argument. Beyond 
4 kpc, this sector represents the Scutum arm, which main- 
tains a positive value of the field, thus implying a field with 
a CC W direction. This r esult is in agreement with the find- 
ings of lHan et al.l ((2006), who reported a CCW field across 
the whole of the Crux-Scutum arm. Our analysis dealt with 
a very limited sample of only 5 pulsars in this sector — one 
of which lies nominally beyond 50 kpc and was therefore 
excluded — and thus the available data did not allow for a 
better probing of the Scutum region. 

Nearly all results from modelling this sector are incom- 
patible with the observed data (see Fig. [6g,h). An excep- 
tion, perhaps, are the field values from the Dipolar-Toroidal 
model, which weakly appear to mirror the magnetic field 
trend — especially for d > 2 kpc. However, the level of 
consistency is low and cannot justify giving any particular 
model the advantage in this case. 
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Figure 7. Image map of the projected (Bi\ ) values for the solar 
neighbourhood (d < 2 kpc). See Fig. [5] for details. 



6.1.4 Norma 

Previous studies of the field direction in the inner Galaxy 
have suggested that the innermost spiral arm, the Norma 
arm, has a CCW field across its whole extent (e.g. Han et 
al. 2002). As it can be seen in the maps of Fig.[5^,c, the field 
inside the arm (5-8 kpc) is mainly negative (apart from a 
very localised positive value around 6 kpc), which implies a 
CCW field. The picture becomes more complex in the Crux- 
Norma interarm region (4-5 kpc), where the field shows an 
abrupt, negative spike in both the plot from the total sample 
and that from only the pulsars below the CP. Above the CP, 
the region is poorly sampled and no safe conclusion can be 
drawn. 

The "spiky" behaviour of the field's profile, which is also 
seen in the Crux sector (Fig.[6ji), implies magnetic fields that 
are an order of magnitude larger than the average value of 
the regular Galactic field (~ 2 ^G). The most likely expla- 
nation for such large values is a high RM gradient between 
pulsars, in conjunction with a small change in the electron 
density (see Eq. I13[l . This is a known phenomenon, where 
groups of pulsars in which the distance between pulsars is 
relatively short display large variations of the RM — and 
sometimes with a change of sign. A likely reason for such 
fluctuations is, again, the turbulent componen t of the field , 
which is at least as strong as the regular field iJonesll 19891 ). 
while another possible reason is the presence of Hu regions, 
which can result in steep electron-density gradients, thus 
affecting the measured RMs . 



6.1.5 Solar Neighbourhood 

The existence of reversals within the solar neighbourhood 
(d < 2 kpc) is a common conclusion of many previous stud- 
ies of the GMF based on pulsar polarisation (see e.g. Lyne 
& Smith 1989, Dickey et al. 2001, Han et al. 2006). Fig. [7] 
shows an image map, similar to that of Fig. [Sfc., focused on 
the solar neighbourhood. At first glance, most longitude sec- 
tors inclu de at least one field re versal between 1 and 2.5 kpc 
distance. iLvne fc Smith! J 19891 ) provided evidence for a re- 
versal between I — 0° and 70°, within d ~ 1 kpc. Their 
find can be confirmed by our image map, which shows at 
least one reversal between d = 1 and 2 kpc in the four, best- 
sampled sectors of that region (see Fig. |4p). F urthermore, 
there is evidence for a similar field pattern in Q4: iFrick et al.l 
(2001) performed a wavelet analysis that suggested a rever- 
sal in Q4, between d — 0.6 and 1 kpc. In that region, the 
only reversal in Fig. [7] is that between I — 285° and 290° . 
However, allowing for the large distance uncertainties borne 
from the uneven and coarse spacing between pulsars, one 
could argue that the rest of the reversals seen in that region 
(up to 2 kpc) support Frick's result. 

6.2 Model Refinement via Parameter 
Optimisation 

In section 16.11 it was shown that nearly all the magnetic- 
field profiles along the selected longitude sectors failed to 
reproduce the field variation derived from actual data. This 
inconsistency could be partly attributed to the chosen values 
of the different parameters that each model depends on. It 
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has been known that parameter optimisation can change 
the model predictions from "being in dramatic disagree- 
ment with the data , to a moderately satisfying agreement" 
|Dickev et al.ll200ll ). 

Therefore, starting from the published parameter values 
for the above four models, we applied a simple fitting routine 
that minimises the reduced-x 2 between the observed and 
model RMs. The latter were calculated via path integration 
along the los of the parallel component of the model mag- 
netic field weighted by the NE2001 density (see Appendix 
for more details). For a set of n RMs and a model with two 
parameters, fei and ki2, we can write: 

2/, i \ 1 P-Mobscrvcd — RM mo del(fc, fe)] 2 , 1rN 

Xn(ki,k 2 ) = -}^ - 2 (15) 

" i-i "observed 

where our analysis considered all 150 RMs derived from the 
observations. 

For multiparametric models, where the number of pa- 
rameters that describe the field are m > 2, the fitting proce- 
dure can be performed via more than two parameters. How- 
ever, the presentation of the Xn {ki, &2 • • • k m ) map poses a 
plotting challenge; hence we decided to attempt to optimise 
only a selected pair of parameters for each model. 

If one takes the published parameters for each model at 
face value, then the above fitting can be done through ad- 
justment factors, e and which simply translate the model- 
specific parameter space to dimensionless multiples of the 
published values: i.e. the published values correspond to 
e = C — !• We chose the fit ranges for e and £ as appro- 
priate for each model and constructed Xn maps that helped 
determine the optimum pair of values for each case. The 
results from such fitting procedure are described in more 
detail below. 

6.2.1 TT model 

In Tinyakov and Tkachev's model, the most obvious mor- 
phological parameters are the pitch angle, p, and the dis- 
tance to the nearest field reversal to the GC, do- Our fit- 
ting approach searched for an optimal value in the ranges 
e : [—0.5,3] and £ : [0,6], which correspond to the 
physical-parameter ranges p = — 8°e : [— 24°, 4°] and do = 
— 0.5£ kpc : [—3, 0] kpc, respectively. The resulting Xn m &P 
of that region is shown in Fig. [8^,. It is clear that there 
is not a unique convergence to a minimum value of x 2 f° r 
any particular pair (p, do). On the contrary, the map reveals 
a complex dependency of the x 2 sphere on the model pa- 
rameters. It is worth noting that a number of local minima 
is observed, but this cannot justify their selection as rep- 
resentative of optimal values, as the former are randomly 
scattered across the map. 

6.2.2 HMR model 

Since the geometry of the HMR model shares many common 
characteristics with that of the TT model, we chose the same 
parameters for the x 2 i- e - the pitch angle, which was 
expressed as p = — 10° e, and the distance to the nearest 
field reversal, do ~ — 0.5023£ kpc. We investigated the Xn 
map corresponding to a parameter space with p : [—30°, 5°] 



and do : [—0.5023, 0] kpc. Fig. [SJa shows a blown-up version 
of that map, focusing on the range where small values of 
X 2 were found. Unfortunately, like in the case of the TT 
model, the resulting map appears complex and without a 
clear convergence. Hence, we could not optimise the HMR 
model's parameters via x 2 -htting between the data and the 
model. 



6.2.3 PS model 

The PS model is clearly a multiparametric one, with a dipo- 
lar and a toroidal component in addition to the disc field. 
Focusing on the field magnitude instead of the morphol- 
ogy, we chose to try and find the pair of values that de- 
scribe best the relative magnitude of the magnetic field in 
the disc compared to the toroidal one. This decision was 
driven by the fact that, although there is a dipolar field in 
this model, its planar component, | B^ lpc ,i e (x,y) |~ 10 -7 - 
10 -8 , is 5-6 times weaker than that of the toroidal field, 
Btoroidni(x,y) |~ 10~ 2 ; and given the fact that most pul- 
sars considered are very close to the GP, with more than 
20% of them with | 6 |< 20°, the -By component, which is 
the quantity of interest, will mostly comprise the spiral and 
toroidal components. 

Our fitting process tried to determine the absolute value 
of the spiral-field's magnitude at the solar radius, ^(r©), and 
the absolute value of the toroidal field maximum, 1.5 kpc 
above the GP, at r : i.e. the value of Br.max (»"©)• For the 
spiral field, we expressed this as fe(r ) = 2e /iG, and for the 
toroidal field, as -Br.max^o) = l-5£ fiG. The best Xn value 
was found in the ranges e : [—0.5, 0] and ( : [8, 12] or, in 
terms of magnetic field magnitude, 6(ro) : [—1,0] /iG and 
Br,max(fo) : [12,18] /iG. Fig. [HJ; shows a contour map of 
that region, where the minimum Xn (~ 3685) corresponds 
to {6(r ),B T ,max(r )} = {-0.35,14.7} fxG. The negative 
sign of the spiral field's magnitude expresses the fact that 
the field direction is opposite to that assumed using the 
published values of Eq. [TOl 

6.2.4 Dipolar-Toroidal model 

In the case of a pure dipolar and toroidal field, the fitting was 
performed on the respective distance-scaled field magnitudes 
expressed by the parameters m and n. Using the published 
values, we explored various ranges of the quantities m = 
245e and n — 4.8£. The x 2 contour map of Fig.[SJi shows the 
region with e : [—0.9, 0.3] and £ : [—0.5, 2.5], where we found 
that the best-fit pair of e and £ corresponds to {m, n} = 
{-36.75 ,izG kpc 3 , 5.37 /iG kpc} ( X l « 3516). 

Table [2] summarises the results from the successful fits of 
the PS and the Dipolar-Toroidal models to the data. Us- 
ing the fitted values, we replotted the model predictions for 
the selected regions discussed in section [6~T1 The results are 
shown in Fig. [9] 

The post-fit plots show a marginal improvement on the 
data-model consistency. In particular, the results from the 
PS model appear to follow the field gradients better, al- 
though this is not the case in terms of the absolute field 
magnitudes. In the region of the Carina arm, for example 
(Fig.[9}I), the post-fit results show better agreement between 
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(b)HMR model 




Figure 8. \n maps of the adjustment factors, t and £, from the parameter-fitting based on (a) the TT model, (b) the HMR model, (c) 
the PS model and (d) the Dipolar— Toroidal model. In all cases, e and £ multiply the published values of the physical parameters of each 
model: i.e. (a) p = — 8°e, the pitch angle, and do = — 0.5f kpc, the distance to the nearest field reversal from the GC; (b) p = — 10°e and 
do ~ — 0.5023f kpc; (c) b(rQ) = 2e fiG, the spiral-field magnitude at the Sun, and -Bt, maxC^©) = l-5£ mGj the maximum toroidal-field 
value, 1.5 kpc above the GP, at tq; (d) m = 244. 7t fiG kpc 3 , the distance-weighted toroidal-field magnitude and n = 4.84f fiG kpc, the 
distance-weighted dipole-field magnitude. 



the model and data points, when compared to the respec- 
tive pre-fit plots. On the whole, between the two models 
for which the fits converged, the PS model gives the most 
consistent picture. 



7 SUMMARY & DISCUSSION 

We used our newly determined RMs to examine particular 
directions in the Galaxy where magnetic field reversal have 
been suggested. We did this by using our measurements to 
determine the radial profiles of the average value of the par- 
allel component of the field along specific longitude sectors: 
the selected directions were those tangential to the Crux, the 
Carina, the Scutum and the Norma arms, and in addition an 
investigation of all directions within 2 kpc of the Sun were 



also examined. Furthermore, we presented an overall picture 
of the magnetic field's direction by means of projecting the 
locations of pulsars that lay close to the Galactic plane and 
plotting a colour-coded image map of the field direction and 
magnitude along 5°-longitude sectors. 

Our work confirmed the existence of a field reversal in 
the Carina-Crux interarm regio n (d = 5-10 k p c), wh ich has 
been derived from the work of iBrown et al. 

I (|2007h . based 

on extragalactic sources and pulsars. Also, a predominantly 
CCW field was derived from both analyses, although our 
work found an additional reversal near 8 kpc, deep into the 
Crux arm. 

Furthermore, we concluded that the Carina arm main- 
tains a CW field with the exception of a negative value at 4 
kpc, which reverses the field to CCW. If this field reversal 
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Figure 9. Scatter plots of the projected as in Fig. [6] The model parameters were adjusted according to the fitting procedure 

described in section [(O] 



is real, it confirms previous claims, bv lHan et all (|2006l ). of 
a "Carina anomaly" somewhere between 3 and 5 kpc. 

Despite the poor sampling of the Galactic region tan- 
gential to the Scutum arm, beyond d = 5 kpc, we confirmed 
the existence of a CCW field, as has been previously re- 
ported. The only pulsar available in the Sagittarius-Scutum 



interarm region suggests that the field direction is retained 
there as CCW, whereas a field reversal, at smaller distances, 
changes the field in the Sagittarius arm to CW. 

Previous studies have indicated that the Norma arm 
may possess a CCW field. If this is indeed true, our findings 
are compatible with such field. However, it is worth noting 
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Table 2. The best parameter values that resulted from fitting 
model RMs from the PS and the Dipolar-Toroidal GMF models 
to n = 150 observed RMs. Each pair of model parameters corre- 
sponds to the minimum Xn found in a scan across the respective 
parameter range. The negative sign in the field magnitude implies 
that the field has an opposite direction to that described by the 
published model parameters. 

PS Dipol.-Toroidal 

Pub. Opt. Pub. Opt. 

6(r ) [ixG] 2 -0.35 

BT.maxO"©) [|UG] 1.5 14.7 

m [fiG kpc 3 ] 245 -36.75 

n [fiG kpc] 4.8 5.37 

X 2 n 3685 3516 



that we derived a CCW direction from only two pulsars with 
known RMs inside the arm itself. 

In many of the scatter plots of the magnetic field against 
distance, there exist strong spikes where the field values are 
many times that of the expected average. This displays the 
vulnerability of this method to samples with neighbouring 
pulsars having completely different RMs. 

Finally, an examination of RMs in the local neighbour- 
hood (d < 2 kpc) revealed a number of reversals in Ql and 
Q4. Although this reinforces the evidence for such reversals 
within 2 kpc, the sparse sampling prevented our analysis 
from determining the locations of the reversals more accu- 
rately; and, hence, it was unable to confirm results from 
previous, more focused studies on the region within a kpc 
from the Sun. 

In an effort to use the new and revised RMs as a test 
bench for large-scale GMF modelling, we implemented a 
simple fitting routine. Using a selected number of exist- 
ing model descriptions and the NE2001 free-electron density 
model, we calculated the RM predicted by these models for 
each of the 150 pulsars by means of integrating the electron- 
density- weighted model magnetic field along the line of sight 
to the pulsar. We found that the overall agreement of model 
RMs and observed data is generally poor. Even our attempts 
to tune key parameters of each model by a x 2 -minimisation 
algorithm were only rewarded by limited success. 

The difficulty to fit even a small sample of well deter- 
mined RMs suggests two possibilties. On one hand, it is 
possible that none of the studied models actually represents 
the true magnetic field structure. On the other hand, small 
scale variations in the ISM and the GMF may alter the ob- 
served RMs sufficiently enough to cause the seen deviations 
from the model expectations. The success of the models to 
explain some of the discussed features suggest that the lat- 
ter is the case. The aim must be therefore, as attempted by 
other authors before (e.g. Brown et al. 2007), to combine 
a sample of RMs as large as possible. This sample should 
consist of both pulsar RMs as well as those of extragalactic 
sources to provide calibration on the largest distances. Even 
then, the information may not be sufficient, and additional 
constraints for fitting the models should be obtained from 
observations other than RM measurements (see e.g. recent 
work by Sun et al. 2007). We will revisit such an approach 
in more detail in a forthcoming publication. 



8 CONCLUSIONS 

We used a novel technique to accurately determine rotation 
measures from 150 pulsars, including 46 new and 104 previ- 
ously published rotation measures. Based on these values, we 
made an effort to investigate — and either confirm or reject 
— previous results on the direction of the regular Galactic 
magnetic field at selected regions. 

We checked four models of the regular Galactic mag- 
netic field for consistency with the data: the Tinyakov- 
Tkachev, the Harari-Mollerach-Roulet, the Prouza-Smfda 
and the Dipolar-Toroidal model. The magnetic-field val- 
ues derived from all the models show large deviations from 
the data. After optimising the Prouza-Srm'da and Dipolar- 
Toroidal models, for which the % 2 fitting algorithm con- 
verged to a unique value, there was a slight improvement on 
the agreement between data and models. The model which 
benefitted mostly from the optimisation was the Prouza- 
Smfda model. However, the magnetic-field-magnitude in- 
consistencies remained at large. A probable explanation for 
these inconsistencies is that the models do not provide the 
full description of the Galactic magnetic field and that a 
more complete model, which will include the short-scale, 
turbulent Galactic- magnetic- field component, is needed. As 
well as lacking a short-scale component, our analysis uncov- 
ered the difficulties intertwined with the effort to fit a global 
magnetic field model to the data. It is almost certain that 
our work will be followed up with more extensive data sets 
and more exhaustive algorithms. However, it remains to be 
seen whether a successful model of the entire Galactic mag- 
netic field is realisable. Perhaps we will have to wait until 
the era of the Square Kilometre Array, with its ability to ob- 
tain rotation measures for a huge number of sources, before 
significant further progress can be made. 
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APPENDIX 



ICordes fc Laziol (120011 ) provide a FORTRAN routine to cal- 
culate the electron density, based on their NE2001 model, 
given a position (x,y,z) in Galactocentric coordinates. The 
DM calulation, when using this routine, is performed by 
integration of the electron density function along the los 
to a particular pulsar. Since the pulsar position, the Sun 
(reference) and the density function are expressed in the 
Galactocentric Cartesian coordinate system, the los has to 
be expressed in the same reference frame. So, the posi- 
tion vector r p to any point along the los to the pulsar is 
parametrised as r p (u) — r© + (r — rQ)u (see Fig. I10[) . 
The integration path along the los is, therefore, expressed 
as r los (u) = r*u = (r - r Q )u. 

Hence, the DM integral is expressed as 



DM 



/ n e (x,y,z) — '" a v ' du 
Jo 



dr\ os (u) 



da 



(A-l) 



and the path integral for the calculation of RM, as 

RM = — I n e (x,y,z)x 
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dx los (u) dy los (u) | a dz los (u) 



B x - 



du 



dt 



da 



(A-2) 



du 



O 
Cl 



> 







s eus 










-.Sun 




a ' 

i - $ 


(t) 


| <v 

o 

Uv %, % 




! ! ! ! Wl'. 


be' 1 * 1 ! ! 1 



■10 



X (kpc) 



Figure 10. Schematic of the position vectors that are used in 
the path integration along each pulsar's los for the calculation of 
the model RMs. The coordinate system used here matches that 
of the NE2001 model: a right-handed Galactocentric Cartesian 
system, originating at the GC, and with the positive y-axis pass- 
ing through the GC and the Sun. The azimuthal angle, 9, in this 
system increases clockwise from the GC-Sun direction. 



where the adopted analytical expression of the GMF vector, 
B(x, y, z), is taken from the models described in section I5TT1 



